New perspective on galaxy clustering as a cosmological probe: general relativistic effects 
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We present a general relativistic description of galaxy clustering in a FLRW universe. The observed redshift 
and position of galaxies are affected by the matter fluctuations and the gravity waves between the source galax- 
ies and the observer, and the volume element constructed by using the observables differs from the physical 
volume occupied by the observed galaxies. Therefore, the observed galaxy fluctuation field contains additional 
contributions arising from the distortion in observable quantities and these include tensor contributions as well 
as numerous scalar contributions. We generalize the linear bias approximation to relate the observed galaxy 
fluctuation field to the underlying matter distribution in a gauge-invariant way. Our full formalism is essential 
for the consistency of theoretical predictions. As our first application, we compute the angular auto correla- 
tion of large-scale structure and its cross correlation with CMB temperature anisotropies. We comment on the 
possibility of detecting primordial gravity waves using galaxy clustering and discuss further applications of our 
formalism. 
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I. INTRODUCTION 

Galaxies are known to trace the underlying matter distri- 
bution on large scales and galaxy redshift measurements can 
therefore provide crucial information about the time evolution 
of large-scale structure of the universe. Over the past decade 
rapid progress has been made in this field, following the ad- 
vent of large galaxy redshift surveys such as the Sloan Digital 
Sky Survey (SDSS; 111]) and the Two Degree Field Galaxy 
Redshift Survey (2dFGRS; H), and much higher precision 
measurements with large survey area have opened a new hori- 
zon for the role of galaxy clustering as a cosmological probe 
(see, e.g., HI). 

However, a critical question naturally arises: is the Newto- 
nian description of galaxy clustering sufficiently accurate on 
large scales, close to the horizon scale at high redshift? Gen- 
eral relativity provides a natural framework for cosmology, 
and the general relativistic description is essential to under- 
stand the formation of CMB anisotropies, since the horizon 
size at the recombination epoch is of order three degrees on 
the sky (see, e.g., ]5|]). Therefore, one would naturally expect 
that a similar relativistic treatment of galaxy clustering needs 
to be considered when using galaxy clustering as a cosmolog- 
ical probe. 

Galaxies are measured by observing photons emitted from 
them and the photon path is distorted by the matter fluctua- 
tions and the gravity waves between the source galaxy and 
the observer. The volume element constructed using the ob- 
served redshift and observed angle is different from the real 
physical volume that the observed galaxies occupy, and the 
observed flux and redshift of the source galaxies are also dif- 
ferent from their intrinsic properties. Therefore, the observed 
galaxy number density is affected by the same perturbations 
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given the total number of observed galaxies, and it contains 
additional contributions from the distortion in the observable 
quantities, compared to the standard description that galaxies 
simply trace the underlying matter distribution 8 g = b S m . 

Furthermore, perturbations are gauge-dependent quantities 
and hence they are not directly observable. For example, 
the matter fluctuation S m computed in the synchronous gauge 
is different from the matter fluctuation S m computed in the 
conformal Newtonian gauge, which diverges on large scales. 
The observable quantities such as observed galaxy clustering 
should be independent of a choice of the gauge condition, 
and this implies that the standard description is incomplete. 
Related to this problem is the linear bias approximation. A 
scale-independent galaxy bias factor b assumed in one gauge 
appears as a scale-dependent galaxy bias factor b(k) in an- 
other gauge. Galaxy formation is a local process and its re- 
lation to the underlying matter density should be well defined 
and gauge-invariant. 

These issues are naturally resolved when we construct the- 
oretical predictions in terms of observable quantities. Here we 
provide a fully general relativistic description of galaxy clus- 
tering in a general Friedmann-Lemaitre-Robertson-Walker 
(FLRW) universe, and as our first application we compute 
the cross correlation of CMB temperature anisotropies with 
large-scale structure. Tracers of large-scale structure contain 
additional contributions from relativistic effects and these ef- 
fects on the cross correlation progressively become significant 
at low angular multipoles at high redshift, since the relativis- 
tic effects are significant at the horizon scale and the horizon 
size decreases with redshift. For a photometric quasar sample 
from the SDSS, we find that the predicted signals are larger 
than the standard method would predict at low angular mul- 
tipoles and its deviation is larger than the estimated cosmic 
variance limit. 

The organization of this paper is as follows. In Sec. Illlwe 
describe our notation for a general FLRW metric and solve the 
geodesic equation for photons. In Sec.|III]we discuss the flue- 
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tuation in luminosity distance that affects the observed flux 
of source galaxies, and we present our main results on the 
general relativistic description of galaxy clustering in Sec.lIVI 
In Sec. [V] we compute the angular correlation of large-scale 
structure and its cross correlation with CMB anisotropics with 
the main emphasis on the systematic errors. Finally, we dis- 
cuss the implication of our new results and conclude with a 
discussion of further applications in Sec. [VI] 



II. GEODESIC EQUATION 

We present our notation for the background metric in an 
inhomogeneous universe and solve the geodesic equation for 
photons emitted from galaxies to derive the relation between 
the source galaxies and the observer. 



A. FLRW Metric 

We assume that the background universe is well described 
by the Friedmann-Lemaitre-Robertson- Walker (FLRW) met- 
ric with a constant spatial curvature, 

ds 2 = g ab dx a dx b = -dt 2 + a 2 (t) g a/3 dx a dx p , (1) 

where ait) is the scale factor and g a p is the metric tensor for 
a three-space. The conformal time t is defined as a dr = dt 
with the speed of light c = 1, and it is related to the comoving 
line-of-sight distance, 

r(r) = a(r )(r - r) = a(r ) / -777 = / tj^t , (2) 

it a(t) J H{z) 

where H = H/a = a/a 2 is the Hubble parameter and the 
dot denotes the derivative with respect to the conformal time. 
The subscript represents that the quantities are computed at 
origin in a homogeneous universe. In a flat universe, the co- 
moving line-of-sight distance is coincident with the comoving 
angular diameter distance. From now on we set a(ro) = 1. 

The metric tensor can be expanded to represent its perturba- 
tions for the spacetime geometry and to describe the departure 
from the homogeneity and isotropy, 

ds 2 = -a 2 (l + 2A)dr 2 ~2a 2 B a drdx a (3) 
+ a 2 [(1 + 2D)g a(3 + 2E a(3 ] dx a dx p . 

We can further decompose the perturbation variables depend- 
ing on their spatial transformation properties as B a = B Q a 
and E a j3 = E Q a p + E^a, where Ef?* is the divergenceless 
tensor. We adopted the convention [6] for the eigenmode Q a 
and Q a p of the Helmholtz equations and assumed there is no 
vector mode. Throughout the paper we use Greek indices to 
represent the 3D spatial components, running from 1 to 3, 
while Latin indices are used to represent the 4D spacetime 
components with being the conformal time component. 

Here we will work with the general representation of the 
metric without fixing gauge conditions (see, e.g., J^, 0, Hi), 



but it often proves convenient to understand our general for- 
mulas in conjunction with other gauges such as the conformal 
Newtonian gauge and the synchronous gauge. The metric in 
the conformal Newtonian gauge (see, e.g., [|9|]) is 

ds 2 = -a 2 (1 + 2tp) dr 2 +a 2 [(1 + 2</>)g a p + 2E% /3 ] dx a dx (3 , 

(4) 

and the metric in the synchronous gauge (see, e.g., lllOlO is 

ds 2 = -a 2 dr 2 + a 2 [g a p + h a/3 ] dx a dx f3 . (5) 

Throughout the paper, we adopt as our fiducial model a flat 
ACDM universe with the matter density Q m = 0.24(£l m ft, 2 = 
0.128), the baryon density tt b = 0.042 (Q, b h 2 = 0.0224), the 
Hubble constant h = 0.73, the spectral index n s = 0.954, 
the optical depth to the last scattering surface s = 0.09, 
and the primordial curvature perturbation amplitude A| = 

2.38 x 10~ 9 at k = 0.05 Mpc" 1 (cr 8 = 0.81), consistent with 
the recent results (e.g., llil fllll ). We use the Boltzmann code 
CMBFast lfl~2tl to obtain the transfer functions of perturbation 
variables. 



B. Temporal component: Sachs- Wolfe effect 

The photon geodesic ir a (A) can be parametrized by an 
affine parameter A, and its propagation direction is then k a — 
dx a /dX, subject to the null equation (ds 2 = k a k a = 0). We 
choose the normalization of the affine parameter, such that 
the time component of the null vector represents the photon 
frequency v, measured by an observer in a homogeneous uni- 
verse [ 13]. The null vector is therefore 

k° = -(l + Sv), k a = --{e a +5e a ), (6) 
a a 

where the unit vector e a is the photon propagation direction 
seen from the observer. The spatial component of the null 
vector is obtained by the null condition and we expanded the 
null vector to the first order in perturbations to represent its 
dimensionless temporal and spatial perturbations 8v and 5e a . 

To the zeroth order in perturbations, the photon frequency 
is redshifted as i> oc 1/a in an expanding universe, and 
the geodesic path is described by d/d\ = (a/P)(d/d\) — 
d T — e a d a = —d/dr. Equivalently the affine parameter \ 
describes the same geodesic path x a (x), but in a conformally 
transformed metric g ab = [p ja) g ab (see, e.g., [ 14] for confor- 
mal transformation). We will put tilde to represent quantities 
in the conformally transformed metric. 

The temporal component of the null vector can be inte- 
grated to obtain the relation between t and \ as 

r-To = x-Xo+ I d X ' 5v( X ') , (7) 
Jxo 

where the subscript o indicates that the affine parameter is 
computed at origin in an inhomogeneous universe. The per- 
turbations of the null vector are related to the metric perturba- 
tions as 

e a 6e a = 8v + A - B a e a - D - E a0 e a e> 3 , (8) 
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by the null equation, and as 

d 



dx 



{5u + 2A) = {A-D)~ (B a{0 + E aP ) e a e 



a o 



(9) 



by the temporal component of the geodesic equation 
(k°-bk b — 0). The vertical bar and the semicolon represent 
the covariant derivatives with respect to g a p and g a b, respec- 
tively. 

Consider a comoving observer of which the rest frame has 
vanishing total three momentum. Its four velocity is u a — 
[(1 — A) /a, v a /a] and the observer measures the redshift pa- 
rameter of a source, 



I = =[-){ l+[5v+A+{v a -B a )eX 



{k a u a ) \a 



(10) 



with the spacetime of the source indicated by the subscript s 
and the bracket representing a difference of the quantities at 
two spacetime points. Using Eq. (0, this relation can be fur- 
ther simplified [15] as 



1 + z. = I — 

a. 



(v a - B a ) e a -A 



(ID 



dr 



(A - D) — (B a \p + E a p) e a e 



where r s = r(z s ) is the comoving line-of-sight distance to 
the source galaxies at z s and v a e a is the line-of-sight peculiar 
velocity. Equation dTTb in the conformal Newtonian gauge is 
known as the Sachs-Wolfe effect 11611 . The first square bracket 
represents the redshift-space distortion by peculiar velocities, 
frame dragging, and gravitational redshift, respectively. The 
first round bracket in the integral also represents the gravita- 
tional redshift, arising from the net difference in gravitational 
potential due to its time evolution for the duration of pho- 
ton propagation, and this effect is referred to as the integrated 
Sachs-Wolfe effect. The last terms in the integral represent the 
tidal effect from the frame dragging and the integrated Sachs- 
Wolfe effect from the time evolution of the primordial gravity 
waves. 

Since the redshift parameter in a homogeneous universe 
is defined as 1/a, we define a quantity Sz that relates the 
observed redshift z s of the source and the redshift of the 
source that would be measured in a homogeneous universe 
as l/a s = (1 + z s )(l — Sz), and note that a a = 1 + fi Sr . 
The redshift l/a s of the source in a homogeneous universe 
is not directly measurable and hence Sz is gauge-dependent. 
One can easily verify that for a coordinate transformation 
t — > t' = t + T, the perturbation in the observed redshift 
transforms as Sz — > Sz' = Sz + TCT, while the observed red- 
shift z s is gauge-invariant. 1 



1 Particular attention needs to be paid to the difference between z and 1/a 
in conjunction with Eq. (TT}. Throughout the paper the redshift parame- 
ter z refers to the "observed" redshift, which is different from the gauge- 
dependent redshift parameter Zh in a homogeneous and isotropic universe, 
defined as 1 + = 1/a . 



C. Spatial components: Gravitational lensing effect 

Metric perturbations, sourced by matter fluctuations and 
gravity waves along the line-of-sight, deflect the photon prop- 
agation direction emitted from galaxies and displace their ob- 
served position on the sky. This effect, known as the gravi- 
tational lensing effect, is described by the spatial components 
of the geodesic equation (k a -,bk b = 0) as 



B L 



2D e a + 2E% e 13 ) 



10 - Sv e a + A\ a - Bp^e 13 - D^ a - E\ 



(12) 



Noting that (d/dx) Sx a = —Se a , the spatial components 
of the geodesic equation can be integrated and expressed in 
spherical coordinates to obtain the angular displacements 



dr 



[(B a - B«) + 2(E a P - Ef)ep\ e e 



r -^)§- e (^-D-B a e a -E a0 e a e0)y (13) 



and 



dr 



[(B a - B°) + 2{E a P - Ef)e ] 



7% smf 



r s - r \ d 



(A-D-B a e a - E aP e a e?) \ . (14) 



Apart from the frame distortion described by I a = {B a — 
B") + 2(E af3 — E" l5 )ep, the gravitational lensing displace- 
ment depends only on the spatial derivative of the metric per- 
turbations, i.e., a constant gravitational potential results in no 
observable effect. 

Since the comoving line-of-sight distance to the source in 
an inhomogeneous universe is (to — r s ) and the source po- 
sition t s is related to the observed redshift z s through a s in 
Eq. ( TTTl l. it can be expressed in terms of r(z) and H(z) in a 
homogeneous universe as 



r = tq - t s = r [z s - (1 



Sz] = 



1 - 



1 



(15) 

where H s = H(z s ). Note that we have expanded the argu- 
ment of r{x) in the square bracket around the observed red- 
shift z s of the source. Finally, using the null equation, the 
radial displacement is then obtained as 



Sr = Xo~Xs + e a Sx a - 



(16) 



St 



dr (A-D-B a e a + E a p e a e ) 



With the full solution of the geodesic equation, the angular 
position § of the source galaxies can be obtained by tracing 
backward the photon path and expressed in terms of observed 
angle n= (0, <j) sin (9) as s= [9 + 86, (0 + 50) sin(0 + S6)}. 
Due to the lensing displacement a unit solid angle |cP§| in 
the source plane is distorted to a unit solid angle |<i 2 n| in the 
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image plane. The amplitude of this distortion is described by 
the convergence k as 



d 2 n 



d 2 s 
and therefore 



/. = ; dr 
/o 



= l-^-5ct>- (cote+ A ] m = I _> ; , . ( 17, 



2 r r. 



esc 6* d^ej I a ) + dg(e e a l a ) + cot g eg, 1° 
V 2 (A - £ - B Q e" - £ Q/3 e Q e' 3 ) 1 , (18) 



where V is the differential operator in two dimensional unit 
sphere. In the literature Eq. (T% is often referred to as the 
gravitational lensing magnification [i. However, the angular 
position s of the source galaxies is not observable; its coordi- 
nate value depends on the choice of gauge condition, while the 
spacetime of the source position is physical. Consequently, 
the convergence k in Eq. (U~8t is gauge-dependent, whereas 
magnification should be a gauge-invariant quantity. In Sec. Hill 
we provide a correct gauge-invariant expression for magnifi- 
cation [i. Note that the gravitational lensing displacements 5r, 
56, and 6<fr are also gauge-dependent. 

The standard Newtonian expression for the convergence 
can be obtained with a few approximations: When the New- 
tonian potential and curvature are constant in time as in an 
Einstein-de Sitter universe, we can replace the total derivative 
d/dr by the partial derivative d r . Integrating by part and ig- 
noring the boundary terms yield the standard form lfl7i [l8ll 
as 



dr 



(r s — r)r 



2r, 



dr 



7 2 -\— 
r 2 dr 

(r s -r)r 



dr 



(19) 



where we have used the Newtonian Poisson equation k 2 <p ~ 
—k 2 ip ~ 4irG a 2 8p. Deep inside the horizon where the 
Newtonian approximation is accurate, there is no gauge am- 
biguity and the gravitational lensing magnification is \i ~ 

\d 2 h/d 2 s\ = 1 + 2 k . 



III. OBSERVED LUMINOSITY DISTANCE 

The observed position and the redshift of source galaxies 
are affected by the matter fluctuations and the gravity waves 
between the source galaxies and the observer, and this relation 
is described by the geodesic equation in Sec. [TT] The observed 
flux of the source galaxies is also affected by the same fluctu- 
ations and this relation is described by the fluctuations in the 
luminosity distance. Here we derive the observed luminos- 
ity distance in an inhomogeneous universe (see 01 Hi for 
earlier derivations). 

Consider a source with intrinsic luminosity L and proper 
radius AR S . The flux T Q and redshift z s of the source are 



measured at origin and the observed luminosity distance is 
defined as 



T> L {z s ) = 



4tt^ 



A S V S 



AR S 



(20) 



where A is the scalar amplitude of the four potential of the 
photons and we have used T oc A 2 v 2 . When the wavelength 
of the photons is shorter than the curvature scale, the prop- 
agation of light rays can be locally described by Maxwell's 
equations, and the governing equations are known as the geo- 
metric optics in curved spacetime (see, e.g., j2lll22Tl ). 

The optical scalar equations are the propagation equations 
of the scalar amplitude 



4- (A a) + - Aad = , 
d\ 2 

and the expansion of the wavevector ■& = k a [a 



1 



dx 



d 2 = ~R ab k a k b 



(21) 



(22) 



where R a b is the Ricci tensor in the conformally transformed 
metric g a i, — (v/a) g a \, . To the zeroth order in perturba- 
tions Eq. d22b has no source term in a flat universe and it 
can be integrated to obtain the expansion of the wavevector 
■Q = 2/(x — Xs — Ax s ), where A% s is related to the size of the 
source. Since the proper radius of the source is AR S = \dt\ 
in a local Lorentz frame, it can be expressed in terms of the 
affine parameter x by considering the photon frequency at the 
source as 



- v s = {k a u a ) s = - 



dt 



a s dx 



(23) 



Note that dt in Eq. d23l is defined in the local Lorentz frame 
of the source. Solving Eq. d2Tb for Aa and using AR S = 
a s v s \Axs\/v s yields the observed luminosity distance as 



T> L {z s 



(1 



)AR t 



Oo 

a* 



exp 



a a (1 + z s ) (xo ~ Xs) 



dx 



1 - 



(24) 



dx 



dr) 



in the limit A% s — » 0, and Sfl is the first order perturbation 
of the expansion of the wavevector that can be obtained by 
expanding Eq. d22l . 

Now to solve for 8i3 we integrate Eq. ( 1221 ) along the zeroth 
order solution 



, Sr) 

dx = 
2 



dr 



(» 



") r 



2r, 



S{R ab k a k b ) 



(25) 



with the source term in the integral 

S(R ab k a k b ) = -k? A- i D - 








B 




i 


T 




E 








\ 




' 3" 





E 

k 2 



5 




e a e' 3 



J a/3 



To the linear order in perturbations the photon geodesic is a 
straight path with small distortion and, with the geodesic path 
x a (x) in Sec. HI] we obtain 



E. 



a/3 



(26) 



Noting that the luminosity distance in a homogeneous uni- 
verse is Dl(z) — (1 + z) r(z) and the comoving line-of-sight 
distance r of the source is related to the affine parameter via 
Eqs. (0 and ( TT5b . the observed luminosity distance 1l(z s ) 
can be written as 1 19] 



Vl{z s ) 
D L {z s ) 



1 + (v a - B a ) s e a - A s 



l + z s 
H s r s 



5z 



dr- 



A 



dr — 

n r s 



{A-D)-{B a]p + E af3 )e a e p 



dr (r ' T)r 5(R ab k a ~k b ) 
2 r„ 



Ti-o 



With the full expression for luminosity distance, the mag- 
nification of a source at observed redshift z is defined as the 
ratio of the observed flux T to the flux of the source that 
would be measured in a homogeneous universe: 



L 



1l 



1 - 2 51 L 



(28) 



We have defined the perturbations in Eq. (f27j as 1l(z) = 
Dl(z)(1 + ST>l), and note that written in terms of observable 
variables 6T>l is gauge-invariant and both 1l(z) and Dl(z) 
are evaluated at the observed redshift z. In the Newtonian 
limit, Eq. d25l ) becomes the convergence k and it is the dom- 
inant factor for 51 l- Therefore, we recover the Newtonian 
expressions 51 l — ~k and /j, ~ 1 + 2 n . 



IV. OBSERVED GALAXY FLUCTUATION FIELD 

Drawing on the formalism developed in Sec. [TT] and [Til] 
we present the expression for the observed galaxy fluctua- 
tion field <5 bs, accounting for all the relativistic effects to 
the linear order. Our formalism is crucial for the theoretical 
consistency and the gauge-invariance of the predictions using 
galaxy clustering as a cosmological probe. To construct the 
observed galaxy overdensity field we start by considering a 
gauge-invariant quantity, the total number N to t of observed 
galaxies. The total number of observed galaxies in a small 
volume described by observed redshift z and observed angle n 
can be formulated in terms of a covariant volume integration 
in a four-dimensional spacetime manifold B23I1 . and it is re- 
lated to the photon geodesic x a (x) via 



-g Up Eabcd U 



dx a dx b dx c 
dz 06 !ty 



dz d6 i 



(29) 



where n p is the physical number density of the source galax- 
ies, the metric determinant is \J—g = a 4 (1 + A + 3D), and 
£abcd = £[abcd\ is the Levi-Civita symbol. 



N tnt = 



r 2 sin 9 
(1 + zfH 



dz dO ' 



l + 3D + v a e a + 2 — 
r 



H ^ T 

r 2 dz 



(30) 



The observed galaxy number density n g is then defined in re- 
lation to the total number of observed galaxies and the ob- 
served volume element as 



r 2 sin 9 



dz d9 i 



5t . (27) 



9 (1 + zfH 

and therefore the observed galaxy number density is 



(31) 



,,, 1 + A + 2D + (v a -B a )e a + E al3 e a e 

- (l + z)^-6z-2^-^5z-5z 

oz Hr 

1 + zdH r n 8r~\ 

- 2k+— — 5z + 2— . (32) 

H dz r 

Given the total number of observed galaxies, the observed 
galaxy number density is affected by the matter fluctuations 
and the gravity waves, since the volume element is con- 
structed in terms of observed redshift and observed angle. 
Equation ( 1291 automatically takes into account the full effects 
of the volume distortion described by the photon geodesic 
equation in Sec. HJ 

However, additional distortions arise due to the intrinsic 
luminosity function dn p /dL of the source galaxies. As de- 
scribed in Sec. Hn] the observed flux of the source galaxies is 
affected by the matter fluctuations and the gravity waves be- 
tween the source galaxies and the observer, and magnification 
of the source galaxy flux changes the observed galaxy num- 
ber density. Given an observational threshold Tthi- in flux at 
origin, the physical number density n p in the above equations 
should be modified as 



■7~tlir 



dTo dLdn^ = [Ahr(1 + 2 ^ (33) 
dj- dL 



where n p (L) is the cumulative (physical) number density of 
the source galaxies brighter than L and Lthr = 4tt£>| (z) Fthr 
is the inferred luminosity threshold for the source galaxy sam- 
ple. For a galaxy sample with dn p j dL oc L~ s , the cumulative 
number density can be expanded as n p (ithr)(l — 5 p 51 l), 
and p — 0.4 (s — 1) is the slope of the luminosity function in 
magnitude. 

Furthermore, since we observe galaxies rather than the un- 
derlying matter distribution, we need to relate the physical 
number density n p of the source galaxies to the matter den- 
sity p m . In the simplest model of galaxy formation, the galaxy 
number density is simply proportional to the underlying mat- 
ter density p m , when p m is above some threshold p t dictated 
by complicated but local process involving atomic physics. 
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z = 3 2 1 0.5 




0.0001 0.001 0.01 

wavenumber k (hMpc -1 ) 

FIG. 1 : Power spectra of perturbation variables computed at z — 
in the conformal Newtonian and the synchronous gauges. Vertical 
lines show the comoving line-of-sight distance (k = l/r(z); light 
gray) at each redshift indicated in the legend and the horizon scale 
(k = Ho ; dark gray) today. Near the horizon scale, even power spec- 
tra of matter fluctuations in two gauges differ dramatically, showing 
that gauge effects are substantial and it is nontrivial to relate perturba- 
tion variables to observable quantities. Two distinct choices of gauge 
conditions cannot be used simultaneously around the horizon scale 
(e.g., Newtonian gauge equations with synchronous gauge transfer 
function outputs from CMBFast or CAMB). 



The matter density at the source galaxy position is related to 
the mean matter density at the observed redshift z as 2 



Pm{x a ) 



Pm(To) 



(1 + 5 m ) = p m {z) [1 + 8„ 



3Sz] 



(34) 

and the mean matter density at the observed redshift is 
p m (z) — (3-ffg /87rG)f2 m (l + z) 3 . The combination (5 m — 
3 Sz) is gauge-invariant and is proportional to the matter den- 
sity at the source galaxy position. Within the linear bias ap- 
proximation, the long wavelength fluctuations of the matter 
density p m at a given point effectively lower the threshold for 
galaxy formation and the galaxy number density can be writ- 
ten as 

n p = n p (z) [1 + 6 (5 m - 3 Sz)] , (35) 
and b is a scale-independent linear bias factor. 3 Equation d35l l 



can be contrasted with the gauge-dependent relation S g — 

b S m . 

Finally, putting all the ingredients together the observed 
galaxy fluctuation field can be written as 



Jobs 



b (S m - 3 Sz) + A + 2D 

x 9 . n 1 + z 
1 + 2 — Sz-2 — — 
oz Hr 

1 + z dH ^ 2 Sr 
H dz r ' 



+ (v a - B a )e a + E^e 13 
Sz — 5z — bp ST>l — 2k 

(36) 



where Sz, Sr, n, and ST>l are given in Eqs. ( fTTT l. (fToT ). ( fl"8l ). 
and d27b . respectively. This equation is the main result of 
our paper. Constructed from the gauge-invariant expressions 
and expressed in terms of observables, this result is gauge- 
invariant. Note that in addition to the scalar contributions, 
Eq. (f36l > includes tensor contributions from the primordial 
gravity waves, mainly from the integrated Sachs-Wolfe effect 
in Sz. 

One remaining ambiguity in computing <5 ODS is the time 
lapse St q at origin, representing the departure from To in a ho- 
mogeneous universe. However, this quantity is independent 
of the position and angle of the source galaxies; In practice 
the mean number density n p (z) of the observed galaxies is 
obtained by averaging n g over observed angle n at a fixed 
observed redshift z and St d is absorbed in the monopole set 
equal n p (z). In the Newtonian limit the dominant contribu- 
tion in Sz is the peculiar velocity V and Eq. d36l > reduces to the 
standard relation for redshift-space distortions 112611 and mag- 
nification bias II27I1 as 



Sstd =bSn 



(37) 



Figure[T]illustrates the theoretical inconsistency in the stan- 
dard method by showing the power spectra of perturbation 
variables computed at z = in the conformal Newtonian 
and the synchronous gauges. The power spectra of matter 
fluctuations in two gauges (solid; synchronous, short dot- 
dashed; conformal Newtonian) noticeably deviate from each 
other well before they reach the horizon scale (dark gray), 
reflecting that theoretical predictions in the standard method 
depend on the choice of gauge conditions. In particular, as we 
observe higher redshift, larger comoving scales (light gray) 
are accessible and the horizon scale is smaller, and therefore 
the systematic errors in the standard methods start to become 
significant on progressively smaller scales. The infrared di- 
vergence shown as the dot-dashed line on large scales is an 
artifact in the conformal Newtonian gauge, while the matter 
fluctuation S m (solid) in the synchronous gauge is also gauge- 
dependent. Theoretical quantities plotted in Fig.[T]are not di- 
rectly observable. 



2 The observed redshift z is related to the expansion parameter a of the 
source galaxy as 1 + z = (1 + 8z)/a . 

3 More general ansatz for Eq. 135} can be obtained by 
generalizing the earlier approach 1 24 l25ll as n p = 
fip(z) exp \b L J yf^g d 4 y (<5 m - 3 Sz)(y) W(x - y)] , where 
VV is a local filter function that cuts off small scale fluctuations, and the 
Lagrangian bias is related to the bias in Eulerian space as b = + 1. 



V. CROSS CORRELATION OF CMB ANISOTROPIES 
WITH LARGE-SCALE STRUCTURE 

As the first application of our formalism, we compute the 
angular correlation of large-scale structure and its cross cor- 
relation with CMB anisotropies. In the standard approach, 
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the observed galaxy fluctuation field is written in the Newto- 
nian limit, and neglecting the additional contributions to the 
observed galaxy fluctuation field results in systematic errors 
in the theoretical predictions. We first introduce the formal- 
ism for computing the angular correlations in Sec. IV Al and 
present the angular auto and cross correlations with the main 
emphasis on the systematic errors in Sec. IV B| 



A. Observed angular fluctuation field 

The observed angular fluctuation field can be obtained by 
integrating Eq. d36l ) along the line-of-sight as 



= I dzP(z)5 ohs (z,n) 



(38) 



with the normalized selection function P(z) of the galaxy 
sample. The selection function P(z) can be obtained by av- 
eraging the observed galaxy number density n g at each ob- 
served redshift slice. Since <5 bs m Eq. d36l l is a linear combi- 
nation of perturbation variables T,- L with different weight func- 
tion Wi(r, n, k), it proves convenient to further decompose 
their functional dependence by 

6 ohs (z, ft) = E / o dr J W ^ ft' k ) T *( k > r ) eik X . 

(39) 

and x = (r, n) in the spherical coordinate. The angular fluc- 
tuation field is often expanded as a function of spherical har- 
monics and the observed angular component is then 

a lm = J d 2 n OA) ^4 (&) = E/ (Sf / dzP ^ 



dr / ^n^(n) W t (r,n) T,(k,r) e lkx .(40) 



For most of the perturbation variables such as 5 m , A, and 
D in Eq. (l36l l. the weight function takes the simple form 
W(r) — 5 D (r — r s ), because they are independent of the 
photon propagation direction and its path. The angular depen- 
dence of the integrand is then carried by the plane wave and 
this functional dependence can be further separated by using 
the partial wave expansion 

e *x = 4% J2i l Ji(kx) Y^(k) F im (n) . (41) 
lm 

The line-of-sight velocity V — v a e a is independent of the 
photon path, but depends on the photon propagation direction; 
the weight function is 



1 d 



W(r) = 8 u (r - r s )(-in ■ k) = -5 u (r - r.) I -— J , 

(42) 

and now it is an operator, acting upon the radial part of the 
plane wave. Note that we have explicitly removed the de- 
pendence on the photon propagation direction by making the 



weight function an operator. The weight function for the weak 
lensing convergence k depends on both the photon path and 
its propagation direction, and it is therefore another operator 
acting upon the angular part of the plane wave: 



W{r) 



2r, r 



' ' M Aa = -ia + l)(Sr4'] - (43) 



2 r s r 



where we have used the relation V 2 Y/ m (n) = — I (I + 

l)Kjm(fi). 

Finally, for the initial conditions described by a Gaussian 
random distribution with A|(fc) oc fc" s_1 , the auto correla- 
tion of large-scale structure can be written as 



Ci 



(a lm ai m ) 



An J^A%k)T^k), (44) 



and the cross correlation of CMB anisotropies with large-scale 
structure is 



.nmb* 



'lm a lm 



An 



dk 



Aj(A)e?(A)2KA), (45) 



where we have defined the angular multipole function of 
large-scale structure in Fourier space 

T l {k)=Y J I dzP{z) f 3 drTi(k,r)Wi(r)ji(kr), 



(46) 

and 9; is the angular multipole function of CMB anisotropies 
(see, e.g., ClUS]). 



B. Angular correlations 

Here we consider a quasar sample without spectroscopic 
redshift measurements used for the cross correlation analysis, 
such as the photometric quasar (QSO) sample [29] obtainable 
from the SDSS. The redshift distribution of the sample is as- 
sumed to have the standard functional form 



P(z) dz oc z a exp 



dz , 



(47) 



with (a, (3, zq)=(3, 13, 3.4). The mean and the peak redshifts 
of the sample are 2.7 and 3, respectively. 

Figure [2] shows the systematic errors in theoretical predic- 
tions of the auto correlation {left) of the QSO sample and its 
cross correlation (right) with CMB temperature anisotropies, 
when the relativistic effects are ignored. Compared to our full 
expression in Eq. d36i l. the theoretical predictions in the stan- 
dard method are computed by using S st d — b S m + (bp — 2)k, 
where S m is the matter fluctuation in the synchronous gauge 
and k is the convergence in the conformal Newtonian gauge. 
We have assumed 6 = 2 and (5p — 2) = 0.1 for the QSO sam- 
ple 1I29I1 . and the full sky coverage of the survey is assumed 
for comparison. 

In the standard approach to modeling the cross correlation 
of CMB anisotropies with large-scale structure, the matter 
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1 10 100 1 10 100 



FIG. 2: Systematic errors in theoretical predictions of the auto cor- 
relation (left) of the QSO sample and its cross correlation (right) 
with CMB temperature anisotropies. Attached bottom panels show 
the mean A^ 2 of the measurements, when only the cosmic vari- 
ance is considered. As in the standard practice, the theoretical pre- 
dictions of the angular correlations are computed by using 8 s td ~ 
b 5 m + (5p — 2) « with 5 m in the synchronous gauge and k in the 
conformal Newtonian gauge, and the angular correlations computed 
with our full expression for <5 bs in Eq. l |36t are compared to the pre- 
dictions with 8 a td- Projection along the line-of-sight suppresses the 
large scale modes, where the matter fluctuations in two gauges in 
Fig. Q] differ substantially. Note that at I = 2 the correct theoretical 
prediction is larger by a factor of 1.8 than that one would incorrectly 
predict in the standard method, and it is 1.2-cr away from the esti- 
mated cosmic variance shown as shaded regions. Since the signal- 
to-noise ratio of the cross correlation measurements is largest at low 
angular multipoles (measurements uncertainties are large at I > 10), 
the systematic errors in the standard method could bias the inferred 
cosmology. 



are underestimated. For the tensor-to-scalar ratio r = 0.1 at 
I = 2, the tensor contribution is ~ 1% of the matter fluctua- 
tion. We emphasize again that <5 bs is gauge-invariant and it 
can be computed in any gauges. 

As opposed to the dramatic contrast seen in Fig.[T]the sys- 
tematic errors in the theoretical predictions seem relatively 
small in Fig. [2] The main reason is the projection effect in 
the angular correlation: each Fourier mode is projected along 
the line-of-sight and the amplitude of C; is largely determined 
by the mode k ~ l/r s , and slightly larger scale mode for the 
cross correlation C ; x due to the cancellation of two spherical 
Bessel functions with different distance scales. This projec- 
tion effect highly suppresses the largest scale modes k ~ l/r s 
of the sample, reducing the dramatic difference in the matter 
fluctuations. For computing the cross correlation with CMB 
temperature anisotropies, one would in practice need to com- 
pute Eq. (f37T > with S m replaced by the combination (6 m — 3 Sz) 
computed in the conformal Newtonian gauge to be consistent 
with the calculation of the convergence k. 

Finally, we comment on the impact of the systematic errors. 
At the lowest angular multipole the accurate theoretical pre- 
diction is about a factor of 1 .8 larger than the standard method 
predicts at the 1.2-cr confidence level with the estimated cos- 
mic variance limit. Since the cross correlation signals decline 
rapidly with angular multipole I, the signal-to-noise ratio of 
the measurements is determined by the estimated cosmic vari- 
ance at I < 10. The lower theoretical predictions in the stan- 
dard method underestimate the cosmic variance, resulting in 
additional A% 2 of a few of the measurements. Considering 
that the current detection significance is at the 3-cr level for 
each galaxy sample [29], these systematic errors could bias 
the inferred cosmology. 



fluctuation in large-scale structure correlates with the inte- 
grated Sachs-Wolfe effect in CMB anisotropies. However, the 
observed fluctuation field 5 bs in Eq. (|36l l contains numerous 
new contributions, including the peculiar velocity, the gravita- 
tional potential, and the integrated Sachs-Wolfe effect, when 
written in the conformal Newtonian gauge. Therefore, when 
computing the cross correlation with the observed galaxy fluc- 
tuation field, the correlations of the new contributions are re- 
quired to be considered in addition to the matter correlation. 
For example, the integrated Sachs-Wolfe effect present in both 
CMB anisotropies and large-scale structure directly correlates 
with each other. 

In the conformal Newtonian gauge, Sz <C S m for the QSO 
sample, since the peculiar velocity, the gravitational poten- 
tial, and the integrated Sachs-Wolfe effect are of the same 
order. Therefore, when <5 b s is computed in the conformal 
Newtonian gauge, the correlation of the matter fluctuation 5 m 
contributes most to C; and C* compared to the other numer- 
ous contributions, and the systematic errors in Fig. [2] arise 
mainly from the difference in 6 m of the conformal Newto- 
nian and the synchronous gauges seen in Fig.Q] However, in 
the synchronous gauge, Sz simply results from the integrated 
Sachs-Wolfe effect due to the absence of the peculiar velocity 
and the gravitational potential, and therefore without account- 
ing for Sz the theoretical predictions in the standard method 



VI. DISCUSSION 

We have developed a fully general relativistic description 
of galaxy clustering as a cosmological probe — we have de- 
rived a covariant expression for the observed galaxy fluctua- 
tion field in a general Friedmann-Lemaitre-Robertson- Walker 
metric without fixing a gauge condition and our formalism 
includes tensor contributions from primordial gravity waves. 
The observed volume element is constructed by using the ob- 
served redshift and observed angle in a homogeneous uni- 
verse, while the real physical volume element given the ob- 
servables needs to be constructed by tracing backward the 
photon geodesic in an inhomogeneous universe. This discrep- 
ancy in the observable quantities results in a significant mod- 
ification of the observed galaxy fluctuation field and provides 
a key clue for understanding gauge issues related to the ob- 
servables. 

As our first application, we have computed the angular auto 
correlation of the photometric QSO sample from the SDSS 
and its cross correlation with CMB anisotropies. The cross 
correlation in the standard method arises from the correlation 
of the integrated Sachs-Wolfe effect in CMB anisotropies and 
the underlying matter fluctuation of the QSO sample. How- 
ever, since there are numerous additional contributions to the 
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observed QSO fluctuation field, the correlations of the addi- 
tional terms with CMB anisotropics need to be considered. 
The dominant contribution to the cross correlation still arises 
from the matter fluctuation and the correct theoretical predic- 
tions are larger at low angular multipoles. The systematic 
errors in theoretical predictions are highly suppressed in the 
angular correlations due to the projection effect, but they can 
result in A\ 2 of a few at low angular multipoles. 

We comment on the possibility to detect primordial gravity 
waves using galaxy samples. Primordial gravity waves affect 
the observed redshift and position of galaxies and therefore 
its effect is imprinted in observed galaxy fluctuation fields. 
We find that the tensor contribution in the cross correlation 
is 1% of the scalar contribution from the matter fluctuation 
at the lowest angular multipoles, when the tensor-to-scalar 
ratio is assumed to be r = 0.1 at I = 2. In general, it is 
extremely difficult to isolate tensor contributions in galaxies, 
because they are completely swamped by scalar contributions. 
One possibility is to cross correlate CMB B-mode polarization 
anisotropics with large-scale structure at low angular multi- 
poles as they should be uncorrected in the absence of primor- 
dial gravity waves. However, it may not be feasible in prac- 
tice, as the parity odd quantities need to be constructed from 
the observed galaxy samples. 

The next application of our formalism is to investigate the 
effect on the three-dimensional power spectrum of galaxy 
samples. Recently, Dalai et al. 113011 showed that the pri- 
mordial non-Gaussianity feature can be probed by the scale- 
dependence of galaxy bias on large scales and this new de- 
velopment has spurred an extensive theoretical and observa- 



tional investigation ll3Tll32l[33ll34l[35ll . However, at this large 
scale, where the primordial non-Gaussianity feature can be 
most sensitively probed, relativistic effects become substan- 
tial and observed quantities are significantly different from 
simple theoretical predictions. Therefore, without proper the- 
oretical modeling of observables, cosmological interpretation 
of these measurements would be significantly biased with the 
current data, and even more so with galaxy samples from fu- 
ture dark energy surveys. Correct modeling of the observed 
power spectrum would not only need to account for the dis- 
crepancy in the observable quantities, but also need to account 
for additional anisotropics arising from the angle dependence 
of the observable quantities (Yoo, Fitzpatrick & Zaldarriaga 
in preparation). 
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